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Abstract 

In this first of three articles on the optical absorption of electrons in half-filled 
Peierls-distorted chains we present analytical results for non-interacting tight- 
binding electrons. We carefully derive explicit expressions for the current 
operator, the dipole transition matrix elements, and the optical absorption 
for electrons with a cosine dispersion relation of band width W and dimer- 
ization parameter 5. New correction ("7y"-)terms to the current operator are 
identified. A broad band-to-band transition is found in the frequency range 
W5 < u < W whose shape is determined by the joint density of states for 
the upper and lower Peierls subbands and the strong momentum dependence 
of the transition matrix elements. 
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I. INTRODUCTION 



Polymers and charge-transfer salts are examples for almost ideal one-dimensional mate- 
rials (Farges 1994). Besides their potential technological applications these systems play an 
important test ground for our theoretical understanding of electron-electron and electron- 
phonon interactions in low dimensions. The theoretical study of the corresponding one- 
dimensional models instead of their three-dimensional counterparts has the further advan- 
tage that much larger systems can be handled numerically, and some of the conceptually 
simplest can even be solved analytically. 

If one wants to show that a model is indeed appropriate for a real material one has 
to calculate experimentally accessible quantities. Analytical calculations often require the 
introduction of further (uncontrolled) approximations while the analysis of data from numer- 
ical approaches is hampered by finite size effects. This makes the assessment of the model 
quality a difficult task. Even if a model can be solved without resorting to uncontrolled 
approximations the interpretation and relevance of the results for real materials remains a 
matter of debate. Theoretical investigations concentrate on ideal systems, e.g., infinitely ex- 
tended and perfectly homogeneous, that cannot be made in practice. Thus it should be kept 
in mind that experimental results may dominantly show effects of disorder and interactions 
with the chain environment which are usually neglected in theoretical approaches. 

Despite these conceptual problems the analysis of basic models remains the first task to 
be accomplished in theory. In principle, models will differ in their experimental predictions 
which will allow to discriminate between them. In addition, model parameters are usually 
fitted by comparison with experiment. 

It is widely believed that polyacetylene and other polymers are basically described by 
independent electrons coupled to a linear chain of atoms (Heeger, Kievelson, Schrieffer and 
Su 1988); (Baeriswyl, Campbell and Mazumdar 1992); (Schott and Nechtschein 1994). For 
a fixed lattice distortion the generic model is the dimerized tight-binding (or Hiickel) model. 
The optical absorption for this independent free electron model was already addressed in 
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the past (Baeriswyl, Harbeke, Kiess, Meier, and Meyer 1983); (Abe 1993). In this work we 
comprehensively derive analytical expressions for the dipole matrix elements and the optical 
absorption in the tight-binding approximation. In contrast to previous works we take special 
care of the derivation of the current and the dipole operator. 

Some charge-transfer salts are understood as strongly correlated one- dimensional electron 
systems for half-filled bands (Mott-Hubbard insulators) (Alcacer, Brau, and Farges 1994). 
The extended Hubbard model for interacting electrons on a distorted chain at half-filling 
is considered appropriate for these materials (Mazumdar and Dixit 1986); (Fritsch and 
Ducasse 1991); (Mila 1995). The theoretical study of these model systems is subject to two 
forthcoming articles (Gebhard, Bott, Scheidler, Thomas, and Koch 1996). 

The paper is organized as follows. In section |I| we set up the Hamiltonian and diagonalize 



it. In section |T| we determine the current operator for tight-binding electrons on a Peierls- 
distorted chain. We show that the distortion gives rise to correction terms that have been 
neglected in the literature. In section |IV| we derive the standard formulae for the optical 
conductivity. In section |V| we present explicit expressions for the optical absorption of free, 
independent tight-binding electrons moving on a Peierls-dimerized chain at half band-filling. 
In section [V| we show how to derive the results via perturbation theory using the regular part 
of the electrical dipole operator. This approach is often preferred in semiconductor optics 
for the perturbational treatment of the external field. A summary closes our presentation. 



II. THE DIMERIZED TIGHT-BINDING HAMILTONIAN 

Free, non-interacting spin- 1/2 electrons in a single band do not show optical absorption 
at finite frequencies because excitations with vanishing momentum necessarily have vanish- 
ing energy. All the optical conductivity is thus concentrated in the Drude weight at zero 
frequency. Hence, interactions must be included to obtain absorption at finite frequencies. 

In highly anisotropic, quasi one-dimensional systems the lattice is unstable against a 
static (Peierls-) distortion at low temperatures since the gain in electron kinetic energy due to 
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the distortion overcomes the lattice deformation energy. We assume that an electron transfer 
can take place between neighboring atoms only (tight-binding approximation). Furthermore, 
the wave vector of the lattice distortion, q, is taken commensurate with the electron filling, 
i.e., q — 2kF — 7r/a where a is the lattice spacing of the undimerized chain. This implies 
that we address the half-filled case where the number of lattice sites L equals the number 
of electrons, N = + = L, n = N/L = 1. 

Under these assumptions the electrons are described by a dimerized tight-binding model, 

l=l,a 

Here, c^^ {ci ^) are the creation (annihilation) operators for an electron of spin a in a Wannier 
orbital centered at site I. < S < 1 describes the effect of the bond-length alternation on 
the electron transfer amphtudes. It imphes t — {ti + 12)/2, and S — {t2 — ^i)/(^i + ^2) > 
if one denotes the hopping amplitudes for the long (short) bond by ti (^2)- Usually, one 
supposes an exponential dependence of the electron transfer amplitudes on the distance r 
between sites, t{r) — ioexp(— (r — a)/X). For ri = a(l +rj/2), r2 — a(l — rj/2) this implies 
^0 = Vtit2 and A = rja/ ln(t2/ti) where r] parametrizes the bond length alternation. Both S 
and 77 are usually small, e.g., S — 0.2 and rj — —0.06 are standard values for polyacetylene 
(Baeriswyl, Campbell, and Mazumdar 1992); (Schott and Nechtschein 1994). Note that 6 
and T] are always opposite in sign. 

As usual the Hamiltonian can be diagonalized in momentum space. We apply periodic 
boundary conditions, and introduce the Fourier transformed electron operators as c^^ = 
\/l7^Ef=i exp{ikla)c'l^^ for the L momenta k = 27rm/(La), m = -{L/2), . . . (L/2) - 1. We 
may thus write 

^(^) = E ^(^) (y^k,a^k,a ~ ^k+n/a,a^k+Tr/a,a) ~ '^^(^) {cl+Tr/a,a^k,a ~ ^k,a'^k+7v/a,a) (2) 
|fc|<7r/(2a),CT 

with the dispersion relation e{k) and hybridization function A(/c) defined as 

e{k) — —2tcos{ka) (3a) 
A{k) ^2tSsm{ka) . (3b) 
4 



The Hamiltonian can easily be diagonalized in fc-space by introducing the new Fermion 
quasi-particle operators a^,^ _|_ for these two bands which are related to the original electron 
operators by 



(4a) 
(4b) 



for \k\ < 7r/(2a). The inverse transformation reads 



(5a) 
(5b) 



where a| + /5| = 1 has to be fulfilled for a canonical transformation. We note the helpful 
relations 



'^k,a'^k,a '^k+n/a,a^k+7T/a,a ~ ('^fc f^k) (^fc,<T,-^fc,o-,- '^k,a, + '^k,cr, 



^k,(7^k,a + ^k+n/a,a^k+'!T/a,a ~ ^k,a,-^k,(7,'- + ^k,a,+^k,, 



^k,a^k+TT/a,a ~^ ^k+ir / a,a^k,(T 



Oil. „ _CLi- 



They allow to write the Hamiltonian in momentum space, eq. (H), as 



m = E K 

\k\<TT/{2a),a 



-^k,a;+ ^k,a;-^k,, 



-e{k){ai-Pt)-A{k)2akPk 
eik)2akPk-Aik)ial-Pl: 



(6a) 
(6b) 

(6c) 
(6d) 



(7) 



We demand that the mixing terms vanish. This gives the following relation between ak 
cos</)fc and Pk = sin0fc 



The mixing amplitudes become 

<^k- Pk~- 



tan(2(/)fc) 



A{k) 



m 

e{k) 



(8) 



m 
m 



1 <k) 



\2\ E{k)^ 



Pk = -^ 



Isgn (A(A:)) 



(9a) 
(9b) 

(9c) 
(9d) 



\2\ E{k)^ 

where sgn(x) is the sign function. Note that ak+n/a = —iPk, Pk+n/a = io-k- This imphes 
that ak+q is complex, in general. 

It is easily seen that the Hamiltonian is diagonal in the new operators. The result is 



f{5) = y: miK^., 



+ ^k,cr,+ ^k,CT-^k 



(10) 



fc|<7r/(2a),CT 

Here, -^E{k) is the dispersion relation for the upper (+) and lower (— ) Peierls band. 



E{k) = dt{kf + ^{kf . 



(11) 



The dispersion of both bands in the reduced zone scheme is shown in figure |I]. The band 
width i^W = At. The model describes a Peierls insulator at half-filling because all states in 
the lower Peierls band are filled, all states in the upper Peierls band are empty, and there 
is the Peierls gap = W5 between the two bands. Optical transitions at finite frequency 
are now possible in the frequency range <Tiu} < W. We set ^ = 1 from now on. 



III. CURRENT OPERATOR AND ELECTRICAL CONDUCTIVITY 



To describe optical transitions we need to include the interaction with an external elec- 
trical field which is only slowly varying in space. Only a frequency dependent field will be 
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considered in the following. In general, a one-dimensional system shall be described by a 
time-independent Hamilton operator. Optical excitations are created by the time-dependent 
perturbation (Maldague 1977); (Mahan 1990) 

Hint^-J-^ (12) 
c 

which couples the external vector potential A{t) to the current operator j of the system (c is 
the speed of light). 

It is not trivial to derive the current operator for a dimerized system. We choose the 
Coulomb gauge and set the scalar potential to zero. Then the electrical field £{t) is given by 
the vector potential ^(t) alone, S{t) = —{l/c){dA{t))/{dt), or, upon Fourier transformation, 

/oo poo fli J 

dte^-'m ; m = / ^e-V(^) , (13) 
-oo J-oo ZTT 

A{u!)/c = £{uj)/{iuj). Both fields are real functions of time, so that they both fulfill f{uj) — 
f*{—uj), f — A,£. The electrical field can be taken constant over lattice distances since 
we are interested in optical excitations in the regime of electron volts. These energies are 
small compared to atomic energies, and the corresponding wave vectors are small compared 
to typical electron momenta. 

The presence of the vector potential changes the amphtudes for the electron transfer 
between neighboring unit cells. In the tight-binding approximation they are given by 



If e \2 

— p+-^ +Ktom(x) 



o VI- . . ^;(x) (14) 

2m \ c ' 

where (— e) is the charge of an electron, and the term in square brackets is the band structure 
Hamiltonian for a single electron. $i(x) are the Wannier functions centered at site We 
can apply a gauge transformation to formally eliminate the vector potential from the band 
structure Hamiltonian. The Wannier functions now acquire a field dependence. 



$i(x) ^ $;(x) 



exp U(e/c) j ' A{x')dx' 



$Kx) (15) 



where is the space coordinate of site I. Since the vector potential is independent of the 
position in the unit cell we may take the field dependence outside the integral for the transfer 
matrix elements. Now that the vector field is along the chain we can write 



ti+i,ii^) = exp {-ie{Ri+i - Ri)A/c) ti+i^i . (16) 

This is the well-known Peierls substitution. 

Note that we tacitly assumed that the Wannier wave functions remain undistorted albeit 
the atoms come closer to or farther away from each other due to the distortion. This implies 
|-R;+i ~ Ri — a\ ^ a as a requirement for the validity of the tight-binding approximation. 
Otherwise the atomic wave functions will be deformed which will ultimately result in a bi- 
atomic, molecular structure. In a dimerized system the distances between neighboring sites 
are given by 

Ri+,-Ri = a{l + v{-iy) ; Ri = a{l-{-iyr]/2) (17) 

where t] is the relative distortion. Formally, we may allow \ri\ < 1 to check our results. On 
physical grounds we have to require |?7| <^ 1 to be consistent with the tight-binding approx- 
imation. The small value t] ^ —0.06 in polyacetylen (Baeriswyl, Campbell, and Mazumdar 
1992); (Schott and Nechtschein 1994) indicates that in real materials the corrections due to 
the 77-terms will dominate those due to Wannier function deformations. The latter correc- 
tions should be exponentially small for small jr^l -values since the wave functions themselves 
decay exponentially as function of distance from the lattice site. 

The Hamiltonian depends on the external field via the field dependent hopping ampli- 
tudes. We expand it in powers of the (weak) external vector potential and obtain 

f (5, A) = m - A + ^T'(5, ^) (7) ' + • • • (18) 

which allows us to identify the operator for the paramagnetic 'particle current^ j, and the 
operator for the diamagnetic field current^ T'{6,t]), as 

l,cr 

f'i6,v) = i-t) E {ct,iA,^ + (1 + (-1)''^) (1 + (-1)'^)' • (19b) 

l,cr 

The operator for the time-dependent perturbation theory in the external field is then iden- 



tified as Hint{t) = —j ■ A{t)/c which is the standard form, equation (IT4 ). The geometrical 



distortion increases (decreases) the electron tunneling between neighboring atoms (5-terms) 
but at the same time decreases (increases) the size of the effective dipole (?7-terms). 

The operator for the total current, jtot, is not only given by j because this quantity is 
not gauge invariant. In fact, the particle current is supplemented by the field current since 
the field itself carries momentum. We write 

]tot= 3 + XA/c (20) 

and demand that this quantity be gauge invariant. This fixes the unknown operator X. 
Applying the same gauge transformation as above we can gauge away the term proportional 
to X, 

itot ^ j{A) . (21) 

We expand eq. ( PPD to first order in A and find from eq. (^) that X = e'^a'^T'{S,T]), as 
expected. 



IV. OPTICAL CONDUCTIVITY 

The dielectric function e{uj) and the coefficient for the linear optical absorption a{uj) are 
given by (Haug and Koch 1990) 

eiu) = 1 + (22a) 

UJ 

aiu) = ^IlM^ (22b) 

where Re{. . .} denotes the real part and rib is the background refractive index. It is sup- 
posed to be frequency independent near a resonance. Hence, the real part of the optical 
conductivity directly gives the absorption spectrum of the system. 

Standard time- dependent perturbation theory gives the Kubo formula for the optical 
conductivity (Maldague 1977); (Mahan 1990). The Fourier transform of the total current 
density and the Fourier transform of the electrical field S{uj) define the optical conductivity 
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Ml. Otot(a;)) xM + Arxe2a(r(5,r7))/L 



(23) 



La 

where (...) means the expectation value in the ground state without the perturbation, and 
is the number of chains per unit area perpendicular to the chain direction. Here, x{'^) 
is the current-current correlation function defined by 



X(w) = -j—^ 
La 



(24) 



where j(t) is the Heisenberg current operator for the unperturbed system. 

The current-current correlation function can be spectrally decomposed in terms of exact 
eigenstates \n) (energy En) of the unperturbed system as 



x(^) 



La 



EKoIjI 



n) 



(25) 



_UJ + {En - Eq) iO- {En - Eq) + i'J _ 

Although 7 = 0"^ is infinitesimal we may introduce 7 > as a phenomenological broadening 
of the resonances at = ±(i?„ — Eq). The real part of the optical conductivity determines 
the absorption spectrum. It is given by 

Im{x(^)} M±7r 



Re{a{uj)} 



LU 



Lauj 



5: |(0|i|n)|' [5 {lu - {En - Eo)) - 5 {lu + {En - Eo))] (26) 



which is positive for all u. 



V. OPTICAL ABSORPTION IN PEIERLS INSULATORS 



According to equation ( PB| ) we have to calculate the transition matrix elements \{n\ 
between the ground state |0) (energy Eq) and all exact excited states \n) (energy En). This 
is easily accomplished for non-interacting electrons. 

The current operator in momentum space reads 



J = -eY^ita (1 + (-1)'5) (1 + (-l)'r^) (c+ ^^^q,, - c+,q+i,. 



I A; I <7r/(3, cr 



r]\ dA{k) 



6 J dk 



(27) 
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It is seen that the current cannot simply be derived from the fc-derivative of the Hamilto- 
nian T{S) in momentum space, eq. (Q). There are additional terms proportional to r] due 
to the fact that the distance between two lattice points is Ar = a(l ± i]/2). This has been 
ignored in previous treatments (Genkin and Mednis 1968); (Cojan, Agrawal, and Flytzanis 
1977); (Baeriswyl, Harbeke, Kiess, Meier, and Meyer 1983); (Abe 1993). 

In terms of the new Fermions the current operator can be split into two terms, one 
which acts within the Peierls subbands while the other induces transitions between them, 

J = iintra + jfntev We obtalu 

ifnter = E ^fnterl^) (^^,.,+0^,^,- + ^^,^,-0^,^,+) (28a) 
|fc|<7r/(2a),cr 

ifntra = E ^fntral^) {K'^,~^k,a,- " afc,a,+afc,a,+) (28b) 
|fc|<7r/(2a),CT 



with the transition matrix elements 



(28c) 



ALa(A:) = ea(l-5^)^^. (28d) 

If we restrict ourselves to half-filling or, more generally, to transitions for positive frequencies 
(interband transitions) we may ignore intraband transitions which only contribute to the 
optical conductivity cr(u;) at = 0. 

Now it is very simple to calculate the optical absorption with the help of the spectral 
decomposition after the Hamiltonian has been diagonalized. One obtains 

Re{aico > 0)} = ^ [Ai„ter(A:)]' 5(cc. - 2E(A;)) (29a) 

|fe|<7r/(2a),<7 

= for W6 < u <W (29b) 

2^' ^(cu2 - {W6y) {W^ - u^) 

where = 4t is the band width, and Afj_ is the number of chains per unit area perpendicular 
to the chain direction. 
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The optical absorption after eq. (^) is shown in figure |^. We plot the dimensionless 
reduced optical conductivity 

ujRe\a(uj > 0)} 

Furthermore we replace the energy conservation 6{uj — 2E{k)) by the smeared function 

^(^) = (31) 

to include effects of phonons and experimental resolution. This is equivalent to the replace- 
ment — + 27 in the formulae for the spectral decomposition of the optical conductivity. 

The optical absorption a^uj) is not only given by the joint density of states for the 
upper and lower Peierls band but also reflects the strong momentum dependence of the 
interband matrix element Aintor(^)- The contribution of the van- Hove singularity a.t u = SW 
is thus enhanced compared to the singularity at u = W. This is the well-known "Umklapp 
enhancement" (Philipps 1966). The effect enforced by the contribution of the ?7-dependent 
terms which have not been taken into account in previous works (Genkin and Mednis 1968); 
(Cojan, Agrawal, and Flytzanis 1977); (Baeriswyl, Harbeke, Kiess, Meier, and Meyer 1983); 
(Abe 1993). 

VI. ELECTRICAL DIPOLE OPERATOR 

In semiconductor physics (Haug and Koch 1990) one often prefers to calculate the optical 
susceptibility x°P*(a;) which is related to the dielectric function via 

e{u) = 1 + Anx'^^iuj) . (32) 

Using eq. (|22aj) and (p3|) we flnd that 



iO'' La „ uj^ 



(33) 

We define the dipole operator P via 
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j = -t P,H (34) 



and insert this definition into eq. (p3D . This allows us to write 

lm{x°^\u;)} = ^ E I (Ol^l^)l' - (En - i?o)) - S {u + (E^ - Eo))] (35) 



n 



which is the result for the optical absorption calculated with the time- dependent perturba- 
tion operator 

Hintit) = -PS{t) . (36) 

In principle we could have used this operator for the time-dependent perturbation theory. 
We can employ eq. (|^) to identify the dipole operator. One finds 

P = _eaj:{l- {-iyv/2) clA<^ = -e ^ RiKA^^ (37) 

l,cj l,a 

which is indeed the result for the dipole operator in second quantization for electrons at 
positions Ri. Unfortunately, this operator is sensitive to boundary conditions and singular in 
the thermodynamical limit. Therefore, it cannot be directly used in any practical calculation. 

It is obvious that the dipole operator still contains a regular contribution for interband 
transitions. To identify this contribution we have to work in momentum space and make 
use of eqs. (PUj ) and ( |28a| ). Then we arrive at 

-Pinter = E f^fntcri^) {j^k,a,+'^k,u- "~ ^fc,cr,-^fc,(7,+) (38a) 
|fc|<7r/(2a),cr 



P /T \ _ 

A'interl'^J ~ 2 



(38b) 



for the interband dipole operator (sgn((5) = — sgn(?7)). Note that it is impossible to give the 
expression for an "intraband dipole operator" because the k-th component of the intraband 
current operator j^tra ^q. ( p8b| ) and the diagonalized Hamilton operator T{6) of eq. (|10|) 
are proportional to each other. In any case such an "intraband dipole operator" would not 
give a contribution to the optical (i.e., finite frequency) absorption, and can thus be ignored 
in the case of a Peierls insulator. 

13 



In the form of equation (^) the interband dipole operator can equivalently well be 
used for the calculation of the optical and higher order susceptibilities for Peierls-distorted 
systems. 

VII. SUMMARY 

In this paper we have presented a detailed analysis of the optical absorption in Peierls 
insulators. We observe the expected signature of a broad band-to-band transition with 
van-Hove singularities at the absorption edges. The strong momentum dependence of the 
transition matrix elements suppresses the singularity at the high-frequency edge compared 
to the one at the low-frequency edge. In contrast to previous investigations we gave explicit 
analytical expressions for the optical absorption and the dipole matrix elements for the tight- 
binding band structure. In addition we found correction ("r7"-)terms to the current operator 
and dipole matrix elements which have been ignored in previous works (Genkin and Mednis 
1968); (Cojan, Agrawal, and Flytzanis 1977); (Baeriswyl, Harbeke, Kiess, Meier, and Meyer 
1983); (Abe 1993). They arise because the distance between two neighboring sites is not 
constant but Ar = a(l ± 7^/2). However, they only slightly influence the final result for the 
optical absorption. 

Unfortunately, a direct comparison to experiment is difficult for two reasons. First, 
disorder effects can inhomogeneously broaden single lines. Hence an experimentally observed 
band can very well be a sign of disorder rather than an argument for a Peierls insulator. 
Secondly, when a residual electron-electron interaction is included in a Peierls insulator (Abe, 
Yu, and Su 1992); (Abe, Schreiber, Su and Yu 1992) one can equally well obtain excitons 
which draw the oscillator strength form the band transitions. Hence, single exciton lines on 
the other hand are not a clear-cut indication against a Peierls insulator either. 
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FIGURES 



FIG. 1. Band structure of a Peierls insulator for S = 0.2. 

FIG. 2. Reduced optical conductivity, aredi^ > 0) in a Peierls insulator for S = 0.2, rj = —0.06. 
A broadening of 7 = 0.011^ has been included. 
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Fig. 1: F. Gebhard et al. 
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Fig. 2: F. Gebhard et aL, 
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